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Abstract. 

We present 2D Monte Carlo radiative transfer simulations of prestellar cores. We consider two types of asymmetry: disk-like 
asymmetry, in which the core is denser towards the equatorial plane than towards the poles; and axial asymmetry, in which the 
core is denser towards the south pole than the north pole. In both cases the degree of asymmetry is characterized by the ratio e 
between the maximum optical depth from the centre of the core to its surface and the minimum optical depth from the centre 
of the core to its surface. We limit our treatment here to mild asymmetries with e = 1.5 and 2.5 . We consider both cores which 
are exposed directly to the interstellar radiation field and cores which are embedded inside molecular clouds. 

The SED of a core is essentially independent of the viewing angle, as long as the core is optically thin. However, the isophotal 
maps depend strongly on the viewing angle. Maps at wavelengths longer than the peak of the SED (e.g. 850 yum) essentially 
trace the column-density. This is because at long wavelengths the emissivity is only weakly dependent on temperature, and the 
range of temperature in a core is small (typically T max /T min Ss 2). Thus, for instance, cores with disk-like asymmetry appear 
elongated when mapped at 850 yum from close to the equatorial plane. However, at wavelengths near the peak of the SED (e.g. 
200 yum), the emissivity is more strongly dependent on the temperature, and therefore, at particular viewing angles, there are 
characteristic features which reflect a more complicated convolution of the density and temperature fields within the core. 

These characteristic features are on scales 1/5 to 1/3 of the overall core size, and so high resolution observations are needed 
to observe them. They are also weaker if the core is embedded in a molecular cloud (because the range of temperature within 
the core is then smaller), and so high sensitivity is needed to detect them. Herschel, to be launched in 2007, will in principle 
provide the necessary resolution and sensitivity at 170 to 250 /im. 
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, 1 . Introduction Bonnor 1956). For example, Barnard 68 has been modelled in 

^"j ■ this way by Alves, Lada & Lada (2001 ). 
/\i • Prestellar cores are condensations in molecular clouds that 

b : are either on the verge of collapse or already collapsing (e.g. However, it is evident from 850 pm continuum maps of 
Myers & Benson 1983; Ward-Thompson, Andre & Kirk 2002). prestellar cores, which essentially trace the column-density 
They represent the initial stage of star formation and their study through a core, that prestellar cores are not usually spheri- 
is important because theoretical models of protostellar collapse call y symmetric (e.g. Motte et al., 1998; Ward-Thompson et 
suggest that the outcome is very sensitive to the initial con- al < 19 "; Kirk et al > 2004 )- Indeed ' statistical analyses of the 
ditions. Prestellar cores have been observed both in isolation projected shapes of a large sample of cores (Jijina, Myers & 
and in protoclusters. Isolated prestellar cores (e.g. L1544, L43 Adams, 1999) suggest that prestellar cores do not even have 
and L63; Ward-Thompson, Motte & Andre 1999) have extents spheroidal symmetry and are better represented by triaxial 
<: 1.5 X 10* AU and masses 0.5 - 35 M Q (see also Andre, ellipsoids (Jones, Basu & Dubinski, 2001; Goodwin, Ward- 
Ward-Thompson, & Barsony 2000). On the other hand, prestel- Thompson & Whitworth, 2002). 

lar cores in protoclusters (e.g. in p Oph and NGC2068/2071) This is not surprising, given the highly turbulent nature 

are generally smaller, with extents ~ 2-4x 10 3 AU and masses f star-forming molecular clouds, and the short time-scale on 

~ 0.05 - 3 M Q (Motte, Andre & Neri 1998, Motte et al. 2001). which star formation occurs (e.g. Elmegreen, 2000). Normally, 

Many authors have modelled prestellar cores with Bonnor- prestellar cores are formed - and then either collapse or dis- 

Ebert (BE) spheres, i.e. equilibrium isothermal spheres in perse - so rapidly that they do not have time to relax towards 

which self-gravity is balanced by gas pressure (Ebert 1955, equilibrium structures. Even in more quiescent environments 

where cores can evolve quasi-statically, the combination of 

Send offprint requests to: D. Stamatellos magnetic and rotational stresses is likely to produce significant 

e-mail: D.Stamatellos@astro.cf.ac.uk departures from spherical symmetry. 
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For instance, SPH simulations of isothermal, turbulent, 
molecular clouds by Ballesteros-Paredes, Klessen & Vazquez- 
Semadeni (2003) show (a) that most of the cores that form are 
transient and non-spherical; and (b) that, despite this fact, the 
column density can, more often than not, be adequately fit- 
ted with the Bonnor-Ebert profile, although the parameters of 
the fit depend on the observer's viewing angle. Likewise, the 
FD simulations of magnetic, isothermal, turbulent, molecular 
clouds reported by Gammie et al. (2003) produce transient, tri- 
axial cores. 

Similarly, evolutionary models of individual prestellar 
cores predict the formation of flattened (oblate spheroidal) 
structures, either due to rotation (e.g. Matsumoto, Hanawa & 
Nakamura 1997), or due to flattening along the lines of a bipo- 
lar magnetic field (e.g. Ciolek & Mouschovias, 1994). Other 
models invoke a toroidal magnetic field to create prolate equi- 
librium cores (e.g. Fiege & Pudritz, 2000). Triaxial structures 
can be generated with a suitable combination of rotation and 
magnetic field. 

It is therefore important to investigate, by means of radia- 
tive transfer modelling, how the intrinsic asymmetries inherent 
in the formation and evolution of a core might translate into 
observable asymmetries on continuum maps of cores. Previous 
continuum radiative transfer modelling of prestellar cores has 
examined non-embedded BE spheres (Evans et al., 2001, 
Young et al., 2003) and embedded BE spheres (Stamatellos & 
Whitworth, 2003a), using ID (spherically-symmetric) codes. 
Zucconi, Walmsley & Galli (2001) have used an approxi- 
mate, semi-analytic method to model non-embedded, magnet- 
ically flattened prestellar cores, in 2D. We have recently de- 
veloped a Monte Carlo code for modelling continuum radia- 
tive transfer in arbitrary geometry, and with arbitrary accu- 
racy. Preliminary results have been presented in Stamatellos & 
Whitworth (2003b,c). 

In this paper we develop continuum radiative transfer mod- 
els of non-spherical cores. Since these models are intended to 
be exploratory, rather than definitive, we consider here only 2D 
models (i.e. we impose azimuthal symmetry so that in spherical 
polar co-ordinates (r, 8, (/>) there is no dependence on (/>). 

Since star formation is a chaotic process, it is not sensible 
to appeal to numerical simulations for the detailed density field 
in a prestellar core. From both observations (e.g. Kirk et al., 
2004) and simulations (e.g. Goodwin et al., 2004), it is clear 
that each core has a unique distribution of gas in its outer en- 
velope, and a unique radiation field incident on its boundary. 
Even if existing simulations are a good representation of the 
real dynamics of star formation, they are not presently able to 
reproduce particular sources, and therefore they can only be 
compared realistically with observations in a statistical sense. 

We will be presenting SEDs and isophotal maps for prestel- 
lar cores formed in SPH simulations of star formation in tur- 
bulent molecular clouds in a subsequent paper (Stamatellos, 
Goodwin & Whitworth, in preparation). However, for interpret- 
ing observations of individual cores, it is more appropriate to 
generate SEDs and isophotal maps using simply parametrized 
models which capture generically the different features we 
might hope to detect. 



In this regard, we have been guided by the observations, 
which indicate that prestellar cores have approximately uni- 
form density in their central regions, and the density then falls 
off in the envelope. If the density in the envelope is fitted with 
a power law, n(r) oc r~ n , then rj ~ 2 - 4 . Here rj ~ 2 is char- 
acteristic of more extended prestellar cores in dispersed star 
formation regions (e.g. L1544, L63 and L43), whereas 77 ~ 4 
is characteristic of more compact cores in protoclusters (e.g. p 
Oph and NGC2068/2071). 

These features are conveniently represented by a Plummer- 
like density profile (Plummer, 1915), 

n{r) = n a - „ ?/2 , (1) 



where no is the density at the centre of the core, and ro is the 
extent of the region in which the density is approximately uni- 
form. 

The Plummer-like density profile is ad hoc, but given the 
transient, non-hydrostatic nature of prestellar cores, and the 
coarseness of the observational constraints, this is unavoid- 
able. It has the advantage of being simple, with only three 
free parameters. Uniquely amongst analytic models, it predicts 
lifetimes, accretion rates, collapse velocity fields, SEDs and 
isophotal maps which agree well with observation (Whitworth 
& Ward-Thompson, 2001; Young et al., 2003). It also repro- 
duces approximately the BE density profile, and the density 
profiles predicted by the ambipolar diffusion models of Ciolek 
& Mouschovias (1994) and Ciolek & Basu (2000). Tafalla et al. 
(2004) use a similar density profile to model the starless cores 
L1498 and L1517B in Taurus-Auriga. 

Furthermore, the Plummer-like profile can be modified eas- 
ily to include azimuthally symmetric departures from spherical 
symmetry. We treat two types of asymmetry. In the first type 
(disk-like asymmetry), we construct flattened cores, using den- 
sity profiles of the form 

■A(i) 2 sin"(0) 



1 +, 



n(r, 6) = no ■ 



O7+D/2 



(2) 



The parameter A determines the equatorial-to-polar optical 
depth ratio e , i.e. the maximum optical depth from the cen- 
tre to the surface of the core (which occurs at 8 = 90°), divided 
by the minimum optical depth from the centre to the surface of 
the core (which occurs at 8 = 0° and 8 = 180°). The parame- 
ter p determines how rapidly the optical depth from the centre 
to the surface rises with increasing 8, i.e. going from the north 
pole at 8 = 0° to the equator at 6 = 90°. 

In the second type (axial asymmetry), we construct cores 
which are denser towards the south pole (8 = 180°) than the 
north pole (8 = 0°), using density profiles of the form 

1 + A (j^f sin" (8/2) 



n(r, 8) - no- 



fo+2)/2 



(3) 



Here the parameter A determines the south-to-north pole opti- 
cal depth ratio e , i.e. the maximum optical depth from the cen- 
tre to the surface of the core (which now occurs at 8 = 180°), 
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divided by the minimum optical depth from the centre to the 
surface of the core (which still occurs at 9 = 0°). The parame- 
ter p again determines how rapidly the optical depth from the 
centre to the surface rises with increasing 9, but now going from 
the north pole at 9 = 0° all the way round to the south pole at 
0=180°. 

In all models the core has a spherical boundary at radius 
/? core . 

For the purpose of this paper, and in order to isolate a 
manageable parameter space, we fix «o = 10 6 crrT 3 , ro = 
2 x 10 3 AU, T] = 2, and R COK = 2 x 10 4 AU. These are typi- 
cal values for isolated cores. We can then explore the effect of 
varying A and p, or equivalently e and p. 

In Section 2, we outline the basic principles underlying our 
Monte Carlo radiative transfer code. In Section 3, we present 
results obtained for cores having disk-like asymmetry; we treat 
both non-embedded cores and cores embedded in molecular 
clouds. In Section 4, we present results for cores having axial 
asymmetry. In Section 5, we summarize our results. 



2. Numerical method & initial system setup 

Our method (Stamatellos & Whitworth 2003a) is similar to that 
developed by Wolf, Henning & Stecklum (1999) and Bjorkman 
& Wood (2001), and is based on the fundamental principle of 
Monte Carlo methods, according to which we can sample a 
physical quantity from a probability distribution using random 
numbers. We represent the radiation field of a source (star or 
background radiation) by a large number of monochromatic lu- 
minosity packets (L-packets). These L-packets are injected into 
the system and interact stochastically with it. If an L-packet is 
absorbed its energy is added to the local cell and raises the lo- 
cal temperature. To ensure radiative equilibrium the L-packet is 
re-emitted immediately with a new frequency chosen from the 
difference between the local cell emissivity before and after 
the absorption of the packet (Bjorkman & Wood 2001). This 
method conserves energy exactly, accounts for the diffuse ra- 
diation field, and its 3 -dimensional nature makes it attractive 
for application to a variety of systems. The code (Phaethon) 
has been thoroughly tested using the thermodynamic equilib- 
rium test (Stamatellos & Whitworth 2003a), and also against 
the benchmark calculations defined by Ivezic et al. (1997) and 
by Bjorkman & Wood (2001). 

The code used here is adapted for the study of cores having 
azimuthal symmetry. The core itself is divided into a number of 
cells by spherical and conical surfaces. The spherical surfaces 
are evenly spaced in radius, and there are typically 50-100 of 
them. The conical surfaces are evenly spaced in polar angle, 
and there are typically 10-20 of them. Hence the core is di- 
vided into 500-2000 cells. The specific number of cells used is 
chosen so that the density and temperature differences between 
adjacent cells are small. 

The L-packets are injected from the outside of the core 
with injection point and injection direction chosen to mimic 
an isotropic radiation field incident on the core (see Fig.[Q. We 
first generate the L-packet injection point on the surface of the 
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L-packet injection point 




Fig. 1. L-packet injection into the core: The packet is injected 
from a random point (9,<p) on the surface of the sphere at 
such an angle (9 m ,<p- m ) as to imitate an isotropic radiation field. 
The frequency of the injected packet is chosen from the Black 
(1994) radiation field. 

core using random numbers < R\,%i e [0,1], 

T = flcore, (4) 

9 = cos _1 (l -2«i) , (5) 
<p = 2n*R 2 ; (6) 

and then the injection direction (Fig.[Q, using random numbers 
^3,^4 e [0,1], 

9 m = cos" 1 (<R\ /2 ) , (7) 
cf> in = 2ifR A . (8) 

Here 9 m is the angle between the normal vector to the tangent 
plane at the point of entry and the packet injection vector, and 
(p m is the polar angle defined on the tangent plane. 

We further assume that the radiation field incident on 
the core is the Black (1994) interstellar radiation field (here- 
after BISRF), which consists of radiation from giant stars and 
dwarfs, thermal emission from dust grains, mid-infrared emis- 
sion from transiently heated small grains, and the cosmic back- 
ground radiation. This is a good approximation to the radiation 
field in the solar neighbourhood but it is not always a good 
choice when studying prestellar cores, since in many cases the 
immediate core environment plays an important role in deter- 
mining the radiation field incident on the core. In this work, 
we simulate the effect of the immediate core environment by 
embedding the core in a molecular cloud which modulates the 
radiation field incident on the core. Another option is to esti- 
mate the effective radiation field incident on an embedded core 
by observing it directly (Andre et al. 2003). 

The dust composition (and therefore the dust opacity) in 
prestellar cores is uncertain, but in such cold and dense condi- 
tions, dust particles are expected to coagulate and accrete ice. 
As in our previous study of prestellar cores (see Stamatellos & 
Whitworth 2003a), we use the Ossenkopf & Henning (1994) 
opacities for a standard MRN interstellar grain mixture (53% 
silicate and 47% graphite) that has coagulated and accreted thin 
ice mantles over a period of 10 5 yr at a density of 10 6 cm -3 . 
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Table 1. Model parameters for cores with disk-like asymmetry 



Model ID 


A 


P 


e 


M corc /M G 


T V (0 = O°) 


1.1 


28 


4 


1.5 


2.0 


94 


1.2 


81 


4 


2.5 


5.1 


94 


1.3 


28 


1 


1.5 


5.0 


94 


1.4 


81 


1 


2.5 


7.3 


94 



e : Equatorial-to-polar optical depth ratio 
M corc : Core mass 

t v (8 = 0°) : Visual optical depth from the centre to the surface of 
the core along the pole (6 = 0°). 



3. Disk-like asymmetry 

3.1. The model 

The density profile in model cores with disk-like asymmetry is 
given by 



6 — -3 



n(r,6) = 10° cm 



1+A (2005au) 2 ^) 

1 + ( >- f " 
1 T V2000AU; 



(9) 



Thus the density is approximately uniform in the centre, and 
falls off as r~ 2 in the outer envelope. The core has a spherical 
boundary at radius R core — 2 x 10 4 AU. The degree of asymme- 
try is determined by A and p, and the values we have treated 
are given in Table 1 , along with 



e - 



ry(6> = 90°) 
MO = 0°) 



(10) 



M core and Ty(9 = 0°). e is the equatorial-to-polar optical depth 
ratio, i.e. the ratio of the maximum optical depth from the cen- 
tre to the surface of the core (which occurs at 9 — 90°) to the 
minimum optical depth from the centre to the surface of the 
core (which occurs at 9 = 0°). Figure 2 displays isodensity 
contours on the x — plane for all the cases in Table 1 . 

3.2. Results: core temperatures, SEDs and images 

In non-embedded cores, the dust temperature drops from 
around 17 K at the edge of the core to 7 K at the centre of the 
core, as previous studies have already indicated (Zucconi et al. 
2001; Evans et al. 2001; Stamatellos & Whitworth 2003a). We 
also find that the dust temperature inside cores with disk-like 
asymmetry is 9 dependent (see Figs. [3}, similar to the results 
of Zucconi et al. (2001). As expected, the equator of the core 
is colder than the poles. The difference in temperature between 
two points having the same distance r from the centre of the 
core but with different polar angles 9, is larger for the more 
asymmetric cores (i.e. those with larger e and/or p; Figs. [3j. 
For example, at half the radius of the core (10 4 AU) the tem- 
perature difference between the point at 9 = 0° (core pole) and 
the point at 9 — 90° (core equator), is 5 - 6 K for the p — 4 
models (Figs. [3^, Wp) but only ~ 2 K for the p — 1 models 
(Figs. |3J;, |5Ji). This temperature difference will affect the ap- 
pearance of the core at wavelengths shorter than or near the 



core peak emission, where the Planck function is strongly (ex- 
ponentially) dependent on the temperature. 

The SED (Fig.|4j of a specific core, for the model parame- 
ters we examine, is the same at any viewing angle, because the 
core is optically thin to the radiation it emits (FIR and longer 
wavelengths). Thus, it is not possible to distinguish between 
flattened and spherical cores, using SED observations, unless 
the core is extremely flattened, so that it is optically thick on the 
equator at FIR and longer wavelengths. Using the Ossenkopf & 
Henning (1994) opacities, the minimum optical depth through 
the model core at 200 pm is 2t\(9 = 0°) x (K200/im/*v) — 0.04, 
so we would need to treat much larger values of e (^ 25) in 
order for the SED to be significantly dependent on viewing an- 
gle. 

In contrast, the isophotal maps of a core do depend 
on the observer's viewing angle. Additionally, they depend 
on the wavelength of observation. Our code calculates im- 
ages at any wavelength, and therefore provides a useful tool 
for direct comparison with observations, e.g. at mid-infrared 
(ISO/ISOCAM), far-infrared (ISO/ISOPHOT) and submm/mm 
(SCUBA, IRAM) wavelengths. We distinguish two wavelength 
regions on which we focus: (i) wavelengths near the peak of the 
core emission (150-250 pm; we choose 200 pm as a represen- 
tative wavelength), and (ii) wavelengths much longer than the 
peak (submm amd mm region; we choose 850 pm as a repre- 
sentative wavelength). In each of the above regions the isopho- 
tal maps have similar general characteristics. 

At 200 pm the core appearance depends both on its tem- 
perature and its column density in the observer's direction. It 
is seen in Figs. 16191 that cores with disk-like asymmetry ap- 
pear spherical when viewed pole-on and flattened when viewed 
edge-on. The outer parts of a core can be more or less luminous 
than the central parts, depending on the core temperature and 
the observer's viewing angle. For example, at close to pole-on 
viewing angles the outer parts of the core are more luminous 
than the inner parts of the core (limb brightening; e.g. Figs. |6j- 
|2l 9 =0°). This happens because the temperature is higher in 
the outer parts and this more than compensates for the lower 
column density (since at wavelengths near the peak of the core 
emission and shorter, the Planck intensity B V (T) depends on 
temperature as B V (T) oc e~"l J , a =const). At other viewing an- 
gles the appearance of the core is determined by a combination 
of temperature and column density effects (Figs. 16191 0=30°, 
60°, 90°). This interplay between core temperature and column 
density along the line of sight results in characteristic features 
on the images of the cores. Such features include (i) the two 
intensity minima at almost symmetric positions relative to the 
centre of the core, on the images at 30°(Fig. f7l9> . and (ii) the 
two intensity maxima, again at symmetric positions relative to 
the centre of the core, on the images at 90°(Figs. I6I7> . (It is 
also worth mentioning that although the characteristic features 
appear in symmetric positions relative to both axes of density- 
symmetry, we should expect deviations from symmetry to arise 
if the radiation field incident on the core is not isotropic.) 

We conclude that isophotal maps at 200 pm contain de- 
tailed information, and sensitive, high resolution observations 
at 200 pm, could be helpful in constraining the core density 
and temperature structure and the orientation of the core with 
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a 



b 



d -> 




Fig. 2. Density distribution on the x — plane (a) for a flat- Fig. 3. Temperature distribution on the x — plane, for the 
tened asymmetric core with equatorial-to-polar optical depth models presented in Fig. |2] We plot iso-temperature contours 
ratio e = 1.5 and p = 4 (model 1.1), (b) for a more flattened from 8 to 18 K, every 2 K. 
asymmetric core, with e = 2.5 and p — 4 (model 1.2), (c) for 
a core with e = 1.5 and p — 1 (model 1.3), and (d) for a more 
flattened core, with e = 2.5 and p — 1 (model 1.4). We plot 
iso-density contours every 10 a5 cm -3 . The central contour cor- 
responds to n = 10 5 5 cm -3 . 
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Fig. 4. SED for the core models in Figs.|5JEl model 1.1 (<? = 
1.5, p - 4; short-dashed line), model 1.2 (e — 2.5, p - 4-; solid 
line), model 1.3 (e = 1.5, p = 1; dash-dot line) and model 1.4 
(<? = 2.5, p - l;long-dashed line). The SED of each core is 
independent of the observer's viewing angle. The dotted line 
on the graph corresponds to the incident/background SED. 

respect to the observer. In Fig. |5] we present a perpendicular 
cut through the centre of the core images shown in Fig. [71 We 
also plot the beam size of the ISOPHOT C-200 camera (90", 
or 9000 AU for a core at 100 pc) and the beam size of the 
upcoming (2007) Herschel (13"or 1300 AU for the 170 /an 
band of PACS; 17"or 1700 AU for the 250 /an band of SPIRE). 
ISOPHOT's resolution is probably not good enough to detect 
the features mentioned above. Indeed, a search in the Kirk et al. 
(2004) sample of ISO/ISOPHOT observations (also see Ward- 
Thompson et. al 2002) does not reveal any cores with such 
distinctive features. However, Herschel should, in principle, be 
able to detect such features in the future. 

In the second wavelength region (submm and mm wave- 
lengths) the core emission is mainly regulated by the column 
density (e.g. at 850 /an, Figs. I10ll3i . Thus, the intensity is 
larger at the centre, where the column density is larger. For the 
same reason the core appears flattened when the observer looks 
at it from any direction other than pole-on. It is also evident that 
the peak intensity of the core is much larger when the core is 
viewed edge-on. Therefore, flattened cores are more prominent 
when viewed edge-on. This introduces a possible observational 
selection effect which should be taken into account when study- 
ing the shape statistics of prestellar cores. Low-mass flattened 
cores are more likely to be detected if they are edge-on, than if 
they are face-on. This is true for optically thin mm and submm 
continuum observations, and also for optically thin molecular 
line observations. For example, when comparing the projected 
shapes of condensations from hydrodynamic simulations with 



Fig. 5. A perpendicular cut through the centre of the core im- 
ages presented in Fig.0for model 1.2 at 200 /an (but also in- 
cluding the background radiation field). In addition we plot the 
beam size of ISO/ISOPHOT (90" = 9000 AU for a core at 100 
pc) and the beam size of the upcoming (2007) Herschel (13" = 
1300 AU for the 170 /an band of PACS and 17"= 1700 AU for 
the 250 /an band of SPIRE). ISO's resolution may not be good 
enough to detect any of the features on the graph, but Herschel 
will have a much better resolution and should be able to detect 
such features in the future. 

the observed shapes using solely optically thin continuum or 
optically thin molecular line observations, one should expect a 
lower number of observed near-spherical cores than indicated 
by the simulations. This may be the reason for the small ex- 
cess of high axis ratio cores in the simulations by Gammie et 
al. (2003) (see their Fig. 9). 

3.3. Embedded prestellar cores 

In the previous section, we studied cores that are directly ex- 
posed to the interstellar radiation field (as approximated by the 
BISRF). However, cores are generally embedded in molecular 
clouds, with visual optical depths ranging from 2-10 (e.g. in 
Taurus) up to 40 (e.g. in p Ophiuchi). The ambient cloud ab- 
sorbs the energetic UV and optical photons and re-emits them 
in the FIR and submm (because the ambient cloud is generally 
cold, r c ] ouc | ~ 20 - 100 K). Therefore, the radiation incident 
on a core that is embedded in a cloud is reduced in the UV 
and optical, and enhanced in the FIR and submm (Mathis et 
al. 1983). Previous radiative transfer calculations of spherical 
cores embedded at the centre of an ambient cloud (Stamatellos 
& Whitworth 2003a), have shown that embedded cores are 
colder (T < 12 K) and that the temperature gradients inside 
these cores are smaller than in non-embedded cores. Andre et 
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Fig. 6. Isophotal maps at 200 //m at viewing angles 0°, 30°, 60° 
and 90°, for a flattened core with equatorial-to-polar optical 
depth ratio e = 1.5 and p — 4 (model 1.1). We plot an isophotal 
contour at 5 MJy sr~' and then from 60 to 110 MJy sr , every 
5 MJy sr" 1 . There are characteristic symmetric features due to 
core temperature and orientation with respect to the observer. 
We note that in these isophotal maps, and in all subsequent 
isophotal maps, the axes (x, y) refer to the plane of sky as seen 



Fig. 7. Same as Fig. [6] but for a more flattened core, with 
equatorial-to-polar optical depth ratio e = 2.5 and p = 4 (model 
1.2). 




Fig. 8. Same as Fig. [6] but for a core with equatorial-to-polar Fig. 9. Same as Fig. [6] but for a core with equatorial-to-polar 
optical depth ratio e — 1 .5 and p — 1 (model 1.3). optical depth ratio e — 2.5 and p = 1 (model 1 .4). 
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Fig. 10. Isophotal maps at 850 /urn at viewing angles 0°, 30°, Fig. 11. Same as Fig. [TO] but for a more flattened core, with 
60° and 90°, for a flattened core with equatorial-to-polar op- equatorial-to-polar optical depth ratio e = 2.5 and p = 4 (model 
tical depth ratio e = 1.5 and p = 4 (model 1.1). We plot an 1.2). 
isophotal contour at 1 MJy sr~' and then from 5 to 50 MJy sr~' , 
every 5 MJy sr~'. The core appears elongated when viewed at 
an angle other than = 0°. 
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Fig. 12. Same as but for a core with equatorial-to-polar Fig. 13. Same as Fig. [TO] but for a core with equatorial-to-polar 

optical depth ratio e — 1 .5 and p — 1 (model 1.3). optical depth ratio e — 2.5 and p = 1 (model 1 .4). 
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al. (2003) also found that the temperatures inside embedded 
cores are lower than in non-embedded cores (assuming that 
they are heated by the same ISRF), using a different approach, 
in which they estimated the effective radiation field incident on 
an embedded core from observations. 



3C I i i i i | i 1 i i I i 
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Fig. 15. Temperature distribution on the x — plane, for the 
same model presented in Figs. 12) and[3j3 (e = 2.5, p — 4, model 
1 .2), but embedded in the centre of an ambient molecular cloud 
with visual extinction Ay = 4. We plot iso-temperature con- 
tours from 8 to 13 K, every 1 K. The embedded core is colder 
than the non-embedded core (Fig.[3p) and the temperature gra- 
dient inside the core is smaller. 



Fig. 14. The effect of the parent cloud on cores. Temperature 
profiles of a non-embedded core (model 1.2; dashed lines), and 
of a core at the centre of an ambient cloud with visual extinc- 
tion Ay = 4 (dotted lines), and Ay = 13 (solid lines). The upper 
curve of each set of lines corresponds to the direction towards 
the pole of the core (0 = 0°), and the bottom curve to the di- 
rection towards the core equator (6 = 90°). The core is colder 
when it resides inside a thicker parent cloud (i.e. when it is illu- 
minated by a radiation field that is weakened at short (< 10pm) 
wavelengths), and the temperature differences between differ- 
ent parts of the core are smaller. Thus, the characteristic fea- 
tures on the isophotal maps at wavelengths near the peak of the 
core SED are weakened. 

Here, we examine the more general case of embedded flat- 
tened cores. We model a core with the same set of parame- 
ters as model 1.2 (p = 4, e = 2.5) but embedded in a uni- 
form density ambient cloud with different visual extinctions 
Ay (Ay = 1.086 ry). The ambient cloud is illuminated by the 
BISRF. In Fig.^] we present the temperature profiles at 9 — 0° 
(core pole; upper curves) and = 90° (core equator; lower 
curves), (i) for a non-embedded core (dashed lines; model 1 .2), 
(ii) for the same core embedded at the centre of an ambient 
cloud with Ay = 4 (dotted lines), and (iii) for the same core 
embedded at the centre of an ambient cloud with Ay = 13 (full 
lines). 

Relative to an non-embedded core, a core embedded in an 
ambient cloud with Ay = 4 is colder and has lower temperature 
gradient (cf. Figs 3b and 15). The isophotal maps are similar 
to those of the non-embedded core (cf. Figs. 7 and 11), but 
the characteristic features are less pronounced. This is because 
the temperature gradient inside the core is smaller when the 



core is embedded (see Fig. ll4l i. For example, at half the radius 
of the core (r = 10 4 AU) the temperature difference between 
the point at 6 = 0° and the point at 6 = 90°, is 5 - 6 K for 
the non-embedded core, but only ~ 1.5 K for the same core 
embedded in an Ay =4 ambient cloud. In Fig.[^]we present 
a perpendicular cut through the centre of the embedded core 
image at viewing angle 30°. It is evident that the features are 
quite weak, but they have the same size as in the non-embedded 
core (Fig. 0, and they may be detectable with Herschel, given 
an estimated rms sensitivity better than ~ 1 - 3 MJy sr~' at 
170-250 yt/m for clouds outside the Galactic plane (dependent 
on cirrus confusion). 

For a core embedded in an ambient cloud with Ay = 13, 
the temperature differences between different parts of the core 
are even smaller (^ 1 K at r — 10 4 AU, see Fig. I14i . but char- 
acteristic features persist (e.g. two symmetric intensity minima 
at 60°and 90°, see FigHTJ. 

Thus, continuum observations near the peak of the core 
emission, can be used to obtain information about the core den- 
sity and temperature structure and orientation, even if the core 
is very embedded (Ay ~ 10 - 20). 

3.4. The effect of a UV-enhanced ISRF on embedded 
cores 

We now examine the effect that a UV-enhanced ISRF has 
on the temperature profiles and isophotal maps of deeply em- 
bedded cores. We consider an ISRF that consists of the BISRF, 
plus an additional component of diluted blackbody radiation 
from stars with = 6000 K or = 10000 K, that illumi- 
nates isotropically the ambient molecular cloud in which the 
core resides. We use a dilution parameter = 10~ 13 , so that 
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the total additional luminosity illuminating the ambient cloud Table 2. Model parameters for cores with axial asymmetry 
is 



(11) 



where <x is the Stefan-Boltzmann constant. This additional flux 
enhances the radiation at A — 0.4 pm incident on the ambient 
cloud by a factor of ~ 3 for = 6000 K, and ~ 30 for = 
10000 K, as compared with the standard BISRF. 

In Fig. ^] we present the temperature profiles at 9 = 0° 
(core pole; upper curves) and 9 = 90° (core equator; lower 
curves) for a core with the same set of parameters as model 
1.2 (p = 4, e = 2.5), embedded in a uniform density ambient 
cloud with visual extinction Ay = 13, illuminated by different 
ISRFs. The bottom pair of curves corresponds to illumination 
by the standard BISRF, and the upper two pairs of curves to il- 
lumination enhanced by a diluted blackbody with T* = 6000 K 
and r* = 10000 K. The core is hotter when the ambient cloud 
is illuminated by a more energetic UV field, by 1 - 2 K for 
the models we examine. The temperature differences between 
different parts of the core seem also to increase, but only by a 
small amount (< 0.5 K). This means that the characteristic fea- 
tures on the isophotal maps at wavelengths near the peak of the 
core emission are not significantly changed. 

For an even more energetic illuminating radiation field 
the enhancement would be larger. For example, the external 
UV radiation field incident on p Ophiuchi is estimated to be 
~ 10-100 times stronger than the BISRF, due to the pres- 
ence of a nearby B2V star (Liseau et al. 1999). In this circum- 
stance, we would expect cores to be hotter and the temperature 
differences may then be sufficiently large (~ 2 K) to produce 
detectable features on isophotal maps at 200 fim. However, if 
the illuminating radiation field is very much stonger than the 
BISRF, it is likely to involve significant contributions from a 
few discrete luminous stars in the immediate vicinity. It will 
therefore be markedly anisotropic and this will produce addi- 
tional asymmetries in the isophotal maps, making their analysis 
more difficult. 

Thus, continuum observations near the peak of the core 
emission may reveal features characteristic of the core struc- 
ture even if the core is more deeply embedded (Ay > 20), pro- 
vided that the radiation field incident on the core is sufficiently 
intense, and provided the effect of discrete local sources can be 
treated. 



4. Axial asymmetry 

4.1. The model 

The density profile in model cores with axial asymmetry is 
given by 



1 +A( 

n(r,ff) = 10 6 cnT 3 — v 



2000 AU 



f sm p (9/2) 



1 + GoOOAu) 



(12) 



Thus - as with model cores having disk-like asymmetry - the 
density is approximately uniform in the centre, and falls off as 
r~ 2 in the outer envelope; the core has a spherical boundary 



Model ID 


A 


P 


e 




T V (0 = 0°) 


2.1 


28 


4 


1.5 


1.4 


94 


2.2 


81 


4 


2.5 


3.4 


94 


2.3 


28 


1 


1.5 


2.4 


94 


2.4 


81 


1 


2.5 


6.3 


94 



e : South-to-north pole optical depth ratio 
Mcore : Core mass 

t v (9 = 0°) : Visual optical depth from the centre to the surface of 
the core along the pole (6 = 0°). 



at radius R cole — 2 x 10 4 AU; and the degree of asymmetry is 
determined by A and p. The values of A and p we have treated 
are given in Table 2, along with 



t v (0 = 180°) 



(13) 



t v (0 = O°) 

Mcore and Ty(6 = 0°). e is now the south-to-north pole optical 
depth ratio, i.e. the ratio of the maximum optical depth from 
the centre to the surface of the core (which occurs at 6 — 180°) 
to the minimum optical depth from the centre to the surface of 
the core (which occurs at 9 = 0°). Isodensity contours on the 
x = plane are plotted on Figure 2 for the model with e = 2.5 
and p — 4. 

4.2. Results: core temperatures, SEDs and images 

In Fig.^5] we present the core density profile on the x — plane 
for model 2.2 (p — 4, e — 2.5), and in Fig.l20lthe corresponding 
temperature profile. The temperature drops from ~ 18 K at the 
edge of the core to ~ 7 K at the centre of the core, and the 
denser "southern" parts of the core are colder. The difference 
between regions of the core with the same r but different 9 is 
larger for the p — 4 models than for the p — 1 models, and 
also larger for the more asymmetric models (<? = 2.5) than for 
the less asymmetric models (<? = 1.5). For example, at half the 
radius of the core (10 4 AU) the temperature difference between 
the point at 9 = 0° (core north pole) and the point at 9 = 180° 
(core south pole) is ~ 7 K for the p = 4, e = 2.5 model, ~ 5 K 
for the p = 4, e = 1.5 model, ~ 4 K for the p = 1, e = 2.5 
model and ~ 3 K for the p = 1, e = 1.5 model. 

As in the case of cores with disk-like asymmetry, these tem- 
perature differences result in characteristic features on isopho- 
tal maps at wavelengths near the peak of the core emission. In 
Fig-|^ we present 200 pm images at different viewing angles. 
The core appears spherically symmetric when viewed pole-on, 
but the effects of axial asymmetry start to show when we look 
at the core from other viewing angles (e.g. 30°, 60° and 90°). 
Comparing with the images at 200 fim for cores with disk-like 
asymmetry (Fig. l6l9> . we see that for cores with axial asymme- 
try there is only one axis of symmetry. Thus, the symmetry of 
the characteristic features is indicative of the underlying core 
density structure. These features contain information about the 
core density, temperature and orientation with respect to the 
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observer, and therefore observations near the peak of the core 
emission are important. The resolution of ISO/ISOPHOT was 
not high enough to detect such features. 

The 850 fim images (Fig. 122b map the column density of 
the core along the line of sight. We point out the similari- 
ties between the maps in Fig. [22] and SCUBA observations of 
L1521F, L1544, L1582A, L1517B, L63 and B133 (Kirk et al. 
2004). Further modelling for each specific core is required to 
make more detailed comparisons. 

As in the case of cores with disk-like asymmetry, the SEDs 
of cores with axial asymmetry are independent of the ob- 
server's viewing angle, because they are optically thin at long 
wavelengths. 

5. Discussion 

We have performed accurate two dimensional continuum ra- 
diative transfer calculations for non-spherical prestellar cores. 
We argue that such non-spherical models are needed because 
observed cores are clearly not spherically symmetric, and are 
not expected to be spherically symmetric. Our models illustrate 
the characteristic features on isophotal maps which can help to 
constrain the intrinsic density and temperature fields within ob- 
served non-spherical cores. They demonstrate the importance 
of observing cores at wavelengths around the peak of the SED 
and with high resolution. 
Our main results are: 

- For the cores treated here, which are optically thin at the 
long wavelengths where most of the emission occurs, the 
SED is essentially the same at any viewing angle. For ex- 
ample, in the case of cores having disk-like asymmetry, the 
SED does not distinguish between cores viewed from dif- 
ferent angles. 

- Isophotal maps at submm wavelengths (e.g. at 850 /mi) are 
essentially column density tracers, whereas maps at far in- 
frared wavelengths around the peak of the core emission 
(e.g. 200 fim) reflect both the column density and the tem- 
perature field along the line of sight. Therefore, sensitive, 
high-resolution observations at 170-250 fim (Herschel) 
combined with long-wavelength observations (e.g. 850 fim 
or 1 .3 mm) can in principle be used to constrain the orien- 
tation of a core and the temperature field within it. 

- If we assume a universal ISRF, then cores embedded in 
ambient molecular clouds are colder than cores directly 
exposed to the ISRF and have lower temperature gradi- 
ents within them. As a result, the characteristic features on 
200 fim isophotal maps are weaker for cores embedded in 
ambient clouds with visual extinction less than Ay ~ 10, 
and may even disappear completely for more deeply em- 
bedded cores (Ay > 20). 

- If the ISRF incident on the ambient molecular cloud is en- 
hanced in the UV region (for example, by nearby luminous 
stars) but still isotropic, then the embedded cores are hotter 
and the temperature gradients inside them may be sufficient 
to produce detectable characteristic features, even in deeply 
embedded cores. 



- The shapes of asymmetric cores depend strongly on the ob- 
server's viewing angle. For example, cores with disk-like 
asymmetry appear more flattened when viewed edge-on. 
Our models also indicate that such cores should be eas- 
ier to detect when viewed from near the equatorial plane. 
This may introduce a selection effect that should be taken 
into account when studying the statistics of the shapes of 
cores, using solely optically thin continuum, or optically 
thin molecular-line, observations. 

- If the ambient ISRF is isotropic, the characteristic features 
on 200 fim maps are symmetric with respect to two axes 
for cores with disk-like asymmetry, and with respect to one 
axis for cores with axial asymmetry. Thus, just the sym- 
metry of these features, could be indicative of the core's 
internal density structure. Lack of symmetry in the features 
could indicate triaxiality, but it could also simply indicate 
that the radiation field incident on the core is anisotropic, 
due to discrete local sources. 

Recently Goncalves, Galli, & Walmsley (2004) have pre- 
sented radiative transfer models of axisymmetric and non- 
axisymmetric toroidal cores, and also models of cores heated 
by an additional external stellar source. In their approach, they 
use a similar but independent Monte Carlo radiative transfer 
code. Their models focus on the effect of different core den- 
sity profiles and anisotropic heating of cores, whereas we focus 
on the effect of observing mildly asymmetric cores at different 
viewing angles. Both studies are in general agreement for sim- 
ilar core models, and are helpful for the study of non-spherical 
realistic prestellar cores. 

We are now extending our study by treating the effect of 
an anisotropic illuminating radiation field. We are also study- 
ing the triaxial molecular cores that result from 3-dimensional 
hydrodynamic simulations. 
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Fig. 16. Same as Fig.|5] i.e. model 1.2, but for a core embedded 
in a uniform molecular ambient cloud with visual extinction 
Ay = 4. The characteristic features visible at different viewing 
angles are weaker than in the case of a non-embedded core. 
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Fig. 17. Same as Fig.|5] i.e. model 1.2, but for a core embedded 
in a uniform molecular ambient cloud with visual extinction 
Ay = 13. There are two characteristic intensity minima visible 
at viewing angles 60° and 90°, but at viewing angles 0° and 30° 
these features are essentially invisible. 
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Fig. 19. Density distribution on the x — plane for a core with 
axial asymmetry, south-to-north pole optical depth ratio e = 
2.5 and p = 4 (model 2.2). We plot iso-density contours every 
10° 5 cirT 3 . The central contour corresponds to n — 10 5 ' 5 cm~ 3 . 



Fig. 18. The effect of an UV-enhanced ISRF on embedded 
cores. Temperature profiles of a core with the same set of pa- 
rameters as model 1.2 (p = 4, e = 2.5), embedded in a uniform 
density ambient cloud with visual extinction Ay = 13, that is 
illuminated by the BISRF (solid lines), by the BISRF plus a 
diluted blackbody of = 6000 K (dotted lines) and by the 
BISRF plus a diluted blackbody of T+ = 10000 K (dashed 
lines). The upper curve of each set of lines corresponds to the 
direction towards the pole of the core (6 = 0°), and the bottom 
curve to the direction towards the core equator (6 = 90°). The 
core is hotter when the illuminating radiation field is enhanced 
at UV wavelengths, but the temperature differences between 
different parts of the core are not significantly increased. 




'• i 

Fig. 20. Temperature distribution on the x = plane, for the 
model presented in Fig.^](e = 2.5, p = 4, model 2.2). We 
plot iso-temperature contours from 8 to 18 K, every 2 K. The 
denser, southern parts of the core are colder. 
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Fig. 21. Isophotal maps at 200 pm at viewing angles 0°, 30°, 
60° and 90°, for a flattened core with south-to-north pole 
optical depth ratio e = 2.5 and p — 4 (model 2.2). We 
plot an isophotal contour at 5 MJy sr~' and then from 60 to 
1 10 MJy sr~', every 5 MJy sr~'. 



Fig. 22. Isophotal maps at 850 fim at viewing angles 0°, 30°, 
60° and 90°, for a flattened core with south-to-north pole op- 
tical depth ratio e = 2.5 and p = 4 (model 2.2). We plot an 
isophotal contour at 1 MJy sr~' and then from 5 to 50 MJy sr~', 
every 5 MJy sr~'. The core appears elongated when viewed at 
an angle other than = 0°. 



